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We develop a formalism to describe the formation of bound states in quantum field theory 
using an exact renormalization group flow equation. As a concrete example we investigate a 
nonrelativistic field theory with instantaneous interaction where the flow equations can be 
solved exactly. However, the formalism is more general and can be applied to relativistic field 
theories, as well. We also discuss expansion schemes that can be used to find approximate 
solutions of the flow equations including the essential momentum dependence. 



I. INTRODUCTION 

The formation of bound states was one of the first problems discussed in quantum mechanics. 
Although the quantum mechanical formalism for nonrelativistic particles with instantaneous in- 
teraction is simple and easy to apply, it is much more difficult to treat the problem in quantum 
field theory as needed for example for relativistic particles. The Bethe-Salpeter [l[ equation can 
be used to sum Ladder diagrams but it is difficult to go beyond [2]. For that reason it is sensible 
to look for alternative approaches. 

In this paper we show how a recently derived exact renormalization group flow equation [3] can 
be used to treat the formation of bound states. As a simple example and to develop the formalism, 
we consider well known problem of nonrelativistic particles with instantaneous interaction. In this 
case we can solve the flow equations exactly and we show that our formalism is equivalent to the 
standard quantum mechanical description although it is based on a fleld theoretic approach. 

In principle, since the flow equation is exact, it can describe also more complicated situations 
for example with relativistic particles, retarded interactions or at nonzero density. However, one 
cannot expect to find exact solutions to the fiow equations, there. The merit of our formalism is to 
provide a convenient starting point for approximations. These do not have to rely on the existence 
of a small coupling constant as in usual perturbation theory and have the potential to capture 
nonperturbative effects. We discuss a generic scheme for the construction of such approximations. 
The resulting formalism is quite general and can be seen as a method to resolve the momentum 
dependence of vertices in renormalization group fiow equations. 

We mention that flow equations for bound states have been investigated previously. A setup 
based on Wegner's flow equation for Hamiltonians jil and the similarity renormalization group JH] 
has been employed to investigate for example a two-dimensional particle with contact potential [6(] . 
Note that the flow equation of Wegner and the similarity renormalization group break space-time 
symmetry explicitly. This is in contrast to Wetterich's formulation of renormalization in terms of 
the effective average action or flowing action F^. 0]. Ellwanger demonstrated that this formulation 
can be used to investigate bound state formation for relativistic theories 0]. By tracing the flow 
of the four-point vertex, Ellwanger calculates for example binding energies for the Wick-Cutkosky 
model and finds good agreement with known results in various limits. As a drawback of his 
approach one might see that properties of the bound states are somewhat hidden in the analytic 
structure of the four-point vertex. This makes is it hard to refine approximations or to investigate 
interactions between bound states or spontaneous symmetry breaking. 

These difficulties can be overcome by introducing auxiliary fields for composite operators as 
will be discussed in more detail below. It is important to adapt the composite fields while solving 
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the flow equation. The reason is that one wants to absorb the essential parts of the four-point 
vertex by means of an Hubbard-Stratonovich transformation at each step in the flow. A first 
proposal, how this can be realized with help of a scale-dependent nonlinear field transformation 
was formulated in ref. 0] and used to study bound state formation in the NJL model. The analysis 
in the present paper is based on a recently derived, simple, but nevertheless exact fiow equation 



[3|]. Although the spirit is similar to the formalism proposed in ref. [9], there are some differences 
in the implementation. 

This paper is organized as follows. In chapter [TI] we introduce our notation and the microscopic 
model for nonrelativistic particles with an instantaneous interaction. In the subsequent section IIIII 
we show how this model can be treated using partial bosonization and an exact flow equation where 
the Hubbard-Stratonovich transformation is kept flxed, i.e. independent of the renormalization 
scale k. In this approach we can integrate the flow equation exactly. However, this solution 
relies on some particular features of the nonrelativistic theory and is difficult to generalize to 
more complicated cases. Therefore we discuss an alternative approach to the problem in section 
IIVI Here we work with a scale-dependent Hubbard-Stratonovich transformation. Again, it is 
possible to integrate the resulting fiow equations exactly in the nonrelativistic case. We argue 
at the end of the section that this formalism is suited much better for the generalization away 
from the nonrelativistic and instantaneous approximation. In section |V] we suggest approximation 
schemes and apply our formalism to a nonrelativistic Yukawa potential. Finally, section IVII gives 
a discussion of a somewhat generalized formalism that can be applied for example in relativistic 
quantum field theory or at nonzero density and we draw some conclusions in section IVIIi 

II. MICROSCOPIC MODEL FOR NONRELATIVISTIC PARTICLES 

Let us consider the problem of two nonrelativistic particles interacting via a potential V{x — y). 
We will use a quantum field theoretic description using the functional integral. The microscopic 
action we employ is 



It is useful to change variables of integration according to x = rjixi + r]2X2, y = xi — X2 with 
rji = Ml/ [Ml + M2), ri2 = M2/{Mi + M2). The interaction term becomes then 




>^ S{ixi)o - {x2)o) 1pl{xi)'lp2{x2)- ( 

We use here = J^^ fx ~ I ^-^o / Coulomb problem the interaction is of the form 



(1) 




(2) 




(3) 



It is also useful to introduce the Fourier transformed fields 





p 



PoXq - px, 
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with 



po 



dpo 



27r 

oo Jp 



(2^ 



(4) 



We also introduce an abbreviation for a composition of fields {ipiip2)- Its position and momentum 
space representation reads 



e-'P-e'<^y{i^iij2){p,q), 



(5) 



p,q 



with {ipiijj2){p, Q) = ipi{rjip + q)ip2{V2P — q)- The conjugate field is 



(^lV'2)*(P,^) 



(6) 
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where the reversed order of tpi and ip2 is convenient when these are fermionic fields. Obviously, 
the first variable (x = (xo,x) or p = {po-,p)) describes the center of mass motion, while the second 
variable {y or q) describes the relative motion of the two particles. The fact that no relative time 
or energy appears is due to the instantaneous approximation for the interaction. For relativistic 
systems this point may be different and the generalization of the formalism to that case will be 
discussed in section IVTl 

Using the momentum representation, the interaction term in the microscopic action becomes 



{'4)i%l)2)*{x,y)V{y) {'4)i'4)2){x,y) 



{i'ii'2)*ip,q)V{q- ^) iipii^2)ip,q' 



In the last line we used 



Jy 



(7) 



(8) 



In this work we investigate the microscopic action in Eq. ([T|) using the functional integral 
formalism. The partition function as a functional of the source field j is given by 



Z[j\ = Di;e 



(9) 



In order to derive an exact fiow equation we modify this expression by introducing an infrared 
cutoff term 



Zk[j] = e 



-iWklj] 



iS^[ip]+iASk flj'ip+c.c.} 



(10) 



where we choose 



A5, riip) 

+i^*2{p) 



2Mi 



Rkif 



Rkif) 



2M2 



i^2{p), 



(11) 
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with Rk — )• oo for /c — )• and — )• for — J- 0. The flowing action is now defined as the modified 
Legendre transform 



{j:iP)Mp) + c.c.} - Wk[j] - A5fc[V], 

Pi=l,2 



(12) 



where ip{p) 
equation 0] 



is now the expectation value. The flowing action satisfies the exact flow 



dkTkiiP] = ^STr(r(2)[v;] +nk)~^dknk. 



(13) 



The operation STr is a trace over both continuous and discrete degrees of freedom such as mo- 
mentum or spin. It also includes a minus sign for fermions. In addition we use in Eq. ()13p the 
matrix of second functional derivatives 



\5^2(_p)/ 



(14) 



and the cutoff matrix 






V 













Rk(f) 

2Mi 





\ 






2M2 



5{p-q), 



(15) 



with 



5{p) = {27:)H{p^)5^^\p). 
For large values of k the flowing action approaches the microscopic action 

hm Tk[^] = S[iljl 

fc— >oo 



(16) 



(17) 



while for A; — )• the flowing action equals the quantum effective action T which is the generating 
functional of one-particle irreducible Feynman diagrams 



hm Tk[M=ni^]- 
fe— s-0 

We note that the flow equation (jlSp can be written as 

dkTM = 4^STYln(r(2)[V;] +7^, 



(18) 



(19) 



with the formal derivative dk that applies to the cutoff term 7^^, only. For reviews of the flow 
equation method see |lCH13l|. 
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III. FIXED HUBBARD-STRATONOVICH TRANSFORMATION 



In this section we use a Hubbard-Stratonovich transformation [14l | to calculate the effective 
bound state propagator. We modify the functional integral over the field ip in Eq. ([9]) by including 
another integral over the (composite) field $ 

Z[j,J] = j Dip e'^'i'^'^^'^'-fp^'^^^^^'^''-^^^''-'''^ 

with 



p,Q,q 



V{q- cf) [^p, ^) - (iPMip, cf)] ■ (21) 

For J = the functional integral over <I> contributes only a multiplicative constant to Z which is 
irrelevant for most purposes. In the combined action S = + 5pb the term quadratic in ip cancels 
and we arrive at 

S = ^ Vl(PO - + ^2(P0 - 

- [ {'^*iP,q)Viq- q')iipi^2)ip,q) + c.c.} 

+ [ '^*{p,q)Viq-q')^p,(f). (22) 
J v.a.a' 



Note that there is a certain ambiguity in the precise from of the Hubbard-Stratonovich transfor- 
mation. Linear transformations of the field $ are always possible. For the Hubbard-Stratonovich 
transformed action in Eq. (|22|) it is now possible to solve some parts of the functional integral. 

Let us consider the functional renormalization group equations for this theory when an infrared 
cutoff of the form in Eq. pip is added for the field ip. For the truncation of the effective action 
we use a vertex expansion. Due to the nonrelativistic dispersion relation and the instantaneous 
interaction, the propagator of the field ip and the Yukawa-like coupling ~ ^*ipi'tp2 + c.c. are inde- 
pendent of the scale k. Also, no additional interaction vertex involving the field tp is generated 
by the renormalization group flow. Indeed, all one-loop Feynman diagrams that contribute to 
the flow equations of the above quantities consist of a closed tour of particles and vanish in the 
nonrelativistic few-body limit. This feature is a manifestation of the fact that the nonrelativistic 
few-body problems decouple in the sence that for example the three-body problem can be solved 
independent of the four-body problem. More general, the flow equations for correlation functions 
that determine the n-body problem are independent of the correlation functions determining the 
n + 1-body problem. For a more detailed discussion of these matters we refer to ref. [l5|. As 
another consequence of the decoupling feature, the flow equation for the propagator of the field $ 
depends only on the propagator of the field ip and the Yukawa-like vertex but is independent of 
the higher vertex functions. This implies that the flow equation for the propagator of ^ can be 
solved exactly. More concrete, we write the term of F^ that is quadratic in $ as 

rf = / <^*iP, q)P^iP, q, q'Mp, <f) (23) 
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with a A;-dependent function P^{p,q,^) denoting the inverse propagator. Its flow equation reads 

dkP^{p,g',q")=dki-i) [ Viq',q)Viq,q") 

Jq 

1 

imp + <i)o - 2m[{vip + 0)"^ + Rkiimp + 0^)] + '■^^^ 

1 



iv2P - q)o - 2m[ivip- 0)"^ + Rkiimp- O)^)] + 

We can write this in a symbohc notation as 

dkP^ip) = -dk{VA-\p)V) (25) 

where the "matrix indices" q etc. and the corresponding integrations are implicit. We use the 
abbreviation 



A^\p,qJ) = i f S('\q-q') 
Jqo 



(26) 



imp + 9)0 - 2m[ivip + q)^ + Rkiimp + of)] + 

1 

iv2P-q)o - 2m[imp- Q)'^ + Rkiimp- o)'^)] + 

= S^'\q- ^) (po - ^[(^iP + 9)' + RkiiViP + q)')] 
1 

- ^^[imp-Q)^ + Rkiimp ■ 

Since the only fe-dependence in this expression comes from R^, the flow equation for P$ can be 
integrated directly. For A — > oo we obtain 

P^Mp) - ^*,a(p) = -VA^'V. (27) 

For P$,A we can use 

P^Aip) = V. (28) 

The propagator 

Gkip) = Pilip) = i-VA-.^V + V)-^ (29) 

describes correlations of two particles with relative momentum to two particles with relative 
momentum q. We use here a matrix notation where q and if are indices and P^^kiP^q^f) is 
inverted as a matrix for fixed value of the center of mass momentum p which is conserved. 

To understand the physical meaning of Gkip) ^® consider field configurations $ for which 
has a pole, or for which 

/ P^,kiP^q,q'Mq')=0 (30) 

Jq' 

holds. In the center of mass frame where p = (po, 0,0,0) and using the reduced mass jj, = 
MiM2/(Mi + M2), this condition can be written as 



/ Viq,^Mq') = 
Jo' 



Po-^if + Rkif)) 



m- (31) 
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For = this is, of course, just the Schrodinger equation for the two-body problem. This 
becomes clear in position space where one obtains 



Po 



1 

2^ 



(-v2 + i?,(-v2))-y(f) 



^>(f) = 0. 



(32) 



In principle, one could obtain the bound state propagator for arbitrary values of the center 
of mass momentum p by inverting the expression for according to Eq. ()29p . To do that one 
would first diagonalize the matrix P$ ^ by finding its eigenvalues Xi{p) and eigenfunctions (pi{p, q). 
This is a problem of similar complexity as solving the Schrodinger equation for the two-particle 
problem. It is clear that for p = (E, 0, 0, 0) and A; = the spectrum of eigenvalues of contains 
as many vanishing eigenvalues as bound states exist with energy E. 

Note, that the on-shell information obtained from solving Schrodingers equation (i.e. the binding 
energies and corresponding eigenfunctions) is not sufficient to determine the bound state propagator 
Gk uniquely. In particular one can add to G^ a term that is regular as a function of the center 
of mass momentum p without changing the poles and therefore the on-shell information. It is 
especially useful to consider the combination 



Gkip) = Gkip) - Ga{p) = Gkip) - V- 



(33) 



This is just the part of Gk that is generated from the fiow equation with the "classical" part 
subtracted. Using Eq. ([29]) this can be written as 



Gk{p) = {A{p)-Vr\ 



(34) 



In this representation it is particularly clear that the poles of Gk correspond to the solutions of 
Schrodingers equation for two particles. 

Moreover, it is clear how to diagonalize this propagator, at least for Rk = 0. To that end we 
note that the inverse propagator has the position space representation 



Gklo{p)=Po 



1 



.p'-i^i-^D-vm 



2{Mi+M2) 2^ 

which is diagonalized, of course, by the solutions of the stationary Schrodinger equation 

1 



Gklo{p)9n{y) = (po 



f - En)9n{y)- 



2(Mi + M2) 



(35) 



(36) 



Here, n is a combined index that labels all quantum numbers. For a spherical symmetric potential 
this includes the radial quantum number as well as those for angular momentum. In the basis 



^(.p,y) = ^4>n(.p)gn{y) 



(37) 



we can write 



iGk=o{p))nm = {PO 



1 



2 (Ml + M2 



P En ) 5nm ■ 



(38) 



For completeness we note the explicit form of the Yukawa coupling term in the effective action in 
this basis 



(39) 



Yl^^ = Y, / {K{p)gmv{mM{p,y) + c.c.] 
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For a constant cutoff function = A;^ the basis of functions Qn diagonalizes also Gk for nonzero 
k. However, the situation becomes more comphcated for other choices of where the basis would 
depend on the center of mass momentum p and the flow parameter k. 

At this point we could in principle undo the Hubbard-Stratonovich transformation and thus go 
back to our original formulation of the theory in terms of the fields ip. This can be done by solving 
the field equation 

§ = » m 

for the field $ as a functional of ip. Plugging this solution $[^] into the action leads to 

rf = rfc[V',$M]. (41) 

In general, the dependence of on <I> is complicated and undoing the Hubbard-Stratonovich 
transformation leads to complicated interaction vertices for the field However, it is quite easy 
to calculate the effective vertex with two incoming and two outgoing '(/'-particles. We note that the 
only contribution to this interaction is given by a tree level diagram involving the propagator Gk- 
The corresponding term in '^^^Iq is 

rit? = - I {M2){p,virgn{yi)V{yi) 

yiy2)gUy2){M2){p,y2) (42) 



(V'iV'2)*(p, y)V{y){ipii^2){p, y)- 



p,y 

This expression has a simple interpretation. The first term is the contribution from bound state ex- 
change processes, while the second term is just the classical term. For an electromagnetic Coulomb 
interaction, the second term can be seen as a contribution from photon exchange processes. 

Let us summarize what we have done in this section. Starting from the microscopic model in 
Eq. ([1]) we performed a Hubbard-Stratonovich transformation and introduced the bilocal field $. 
In this "partially bosonized" language it was possible to find a solution of the fiow equation for the 
propagator of the field $ in a closed form. Note however, that this solution relies on some features 
particular to the nonrelativistic few-body problem. First, the description of the interaction between 
particles in terms of a instantaneous interaction potential V{xi — X2) is usually not possible for 
relativistic problems. In relativistic quantum field theory, interactions are mediated by exchange 
particles such as the photon. In momentum space the resulting interaction term has nontrivial 
frequency- and momentum dependence and is subject to renormalization group modifications for 
example due to a "running" fine-structure constant a or charge e. These features make it difficult 
to absorb the interaction term by an appropriate Hubbard-Stratonovich transformation. Even if 
this was possible, solving the flow equations would be more difficult than in the nonrelativistic case 
presented above. Both the propagator for the field ip and the Yukawa interaction ~ ^*'ipi^2 + c.c. 
might have non-vanishing fiow equations. These equations as well as the fiow equations for the 
propagator of $ would depend on higher-order vertex functions such that no solution in a closed 
form can be expected. 

Although one cannot expect to find analytic solutions by transferring the above calculation 
to a relativistic field theory, one might ask whether it can be helpful for finding approximate 
solutions. One could make a truncation of the flowing action F^ for example in terms of a derivative 
expansion and consider the flow equations in the theory-subspace spanned by the flnite number 
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of operators included in this truncation. Indeed, we will argue below that such truncations can 
lead to good approximate solutions. The formulation in the present section has one drawback 
for approximate solutions, however. As already discussed above, interactions in relatvistic field 
theories are usually mediated by exchange fields such as the photon. The renormalization of 
the corresponding couplings and propagators can be calculated most efficient in a formulation 
of the theory which directly takes the exchange fields as propagating fields into account. Other 
formulations where these fields (as for example the photon) are "integrated out" might be equivalent 
in principle, but are usually much harder to treat by approximate methods. The reason is that the 
essential momentum- and frequency dependence is often hidden in such formulations. 

Transfered to the example of a nonrelativistic theory it would be useful to have a formulation 
where the two terms in Eq. ()42p are treated in different ways. While the term in the first line 
(the contribution from bound state exchange processes) is treated most efficient in terms of the 
Hubbard-Stratonovich field this is different for the second term. For a relativistic field theory 
this corresponds to the contribution from photon or other exchange processes and it is therefore 
most efficient to write it in terms of this exchange field. In other words, for the nonrelativistic 
theory, the Hubbard-Stratonovich transformation should be constructed such that the effective 
action Tk=o contains as an explicit contribution only the classical term (the second term in Eq. 
([12]) ). All additional terms should be described by bound state exchange processes. For the flowing 
action one can use the same prescription. To that end it is necessary to work with a scale-dependent 
Hubbard-Stratonovich transformation, however. An exact flow equation that can be used for this 
purpose was derived m i and wih be discussed m the next section. 



IV. SCALE DEPENDENT HUBBARD-STRATONOVICH TRANSFORMATION 

In this section we investigate the bound state problem using a ^-dependent version of the 
Hubbard Stratonovich transformation. Instead of Spi, in Eq. ([2T]) we employ 



5pb= / [^*{p)-{i^M*(j^)VQl\p)] Qip) 

Jv 



'V (43) 

[^(p)-Qa'(p)^(V'iV'2)(p)] . 

Again we use a matrix notation where the summation over the index q etc. is left implicit. In Eq. 
([I3|) the matrix Qip) is ^-dependent, while Q\{p) equals Q{p) for A: = A but is independent of k. 
We choose 



lim QA{p,q,q') = oo 6^^\q- q'), 
A— s-oo 

lim Ql^{p,q,q') = 0. 

A^oo 

In this limit the combined action 5 = 5^ + S'pb becomes 

S = I i^lipo- TTTr^^ ) ^1 + V'2 (po - Tirrf ] ^: 



(44) 



p 



2Mi" / ' -^M"" 2M2" ' 



{M2Y{p,q)V{q- q)i'4'iil^2)ip,q) 



{^*{p,q)ViQ- 'f){i^M{p,<f) + c.c.} 



p^m 



+ / ^*{p,q)Q{P,q.(imP,^)- (45) 
Jp4,q' 
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The field $ is very "massive" for large A. The action in Eq. (j45p is the starting point for the 
flow of the functional Tk[il^, ^] for k = A. In ref. 0] we showed that the flowing action for a scale 
dependent Hubbard-Stratonovich transformation of the form in Eq. (j43p satisfies the exact flow 
equation 

dkVk = ^STr(rf +7^fc)-l(afc7^,.-7^fc(5fcQ-l)7^fe) 

-l^'i\dkQ'')r^'> + (46) 

where 7^ is a field independent constant that is irrelevant for most purposes and will be dropped 
from here on. The flow equation (j46p holds for fixed (/c-independent) field ^. For our purpose it is 
useful to perform a fc-dependent linear transformation on the field 

^*ip, q) = ^"{p, ^) (e^^^(P)) (g', 0, (47) 



with an exponentiated matrix e^^''^^ with "indices" q,(f. This gives the flow equation 
^fcrfcli, = dkTk\^ - [ -Ip^dk^{p,q)\ 



,q5<^{p,q) 

{d,^*ip,q))\^T^^ 
p,g '^6^*{p,q) 

-§i(a,M,)$-c|.*(a,Mt)^. (48) 
In the last equation we used an obvious matrix notation and the abbreviations 



k 



(49) 



and functional derivatives are now taken with respect to It is important to note that the 
matrices 9fcQ, d^Q and have entries only in the <I>-<I>-block. For example, Eq. (j49p implies 
a transformation of the cutoff term for the composite bosons but not for the fundamental fields 
ipi,ip2- For simplicity we drop the hats at most places below. 

To investigate the implications of the flow equation (|48p we use again a truncation of the 
flowing action in terms of a vertex expansion. Due to the nonrelativistic dispersion relation and 
the instantaneous interaction, this leads to exact flow equations for the considered n-point functions 
since their flow equations decouple from higher vertex functions. More concrete, we choose as our 
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truncation 



(V'iV'2)*(?', q)>^^{p, q, g )(V'iV'2)(p, g') 



+ / <^*{p,q)P^{p,q,q'MP,<f)- (50) 

Once again, the propagators for the fields ipi,ip2 remain unmodified by the renormahzation fiow, 
while the functions A^, and P$ are /c-dependent objects. For k = A they have the initial values 

Av>,A(p,g,g') = 0, 
h<s>,A{P^Q^Q) = y{q-^)^ 

P'S'AiP^ Q) = Qdq. q') = oo [q- q'). (51) 

Projecting the flow equation ()48p onto the truncation in Eq. ()50p we find for A^, = the flow 
equations (again in matrix notation and suppressing the argument p at several places) 

dk\^ = dk{VA-'V) + h^idkQ-')h^, 
dkh^ = dkih^A^^V) - P^idkQ-^)h^ 

~{dkMl)h^, 

dkP^ = -dk{h^A-^h^) - P^{dkQ-^)P^ 

~{dkMl)P^ - P^idkMk). (52) 

As discussed in ref. 0] we can now use our freedom in the choice of dkQ~^to enforce d^Xip = 0, i.e. 

X^{p,q,q)=0 for all A;. (53) 

This fixes dkQ~^ = —h'^^V{dkA'^^)Vh^^. In addition, the A;-dependent field rescaling determined 
by the matrix Mk{p) can be chosen arbitrary as well such that we can use it to enforce dkh^{p) = 0, 
i.e. 



h^iPi q) = ^{q — q) for all k. (54) 



In summary, we use 



dkQ-^ = -dkA-\ 
dkMk = {dkAl^)\P^ + V)\ 

dkMl = {p^ + V){dkA^^). (55) 

For P$ this leads to the flow equation 

dkP^ = -V{dkA-^)V + P^{dkA^^)P^ 

-(P$ + V){dkA-')P^ - P^{dkA^'){P^ + V) 
= -{P,i. + V){dkA-'){P,i. + V) (56) 
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which is solved by 

P$ = Ak - V. (57) 

Using the definition of Aj. in Eq. ([2^]) . we see that, as expected, the zero crossings of P$ at /c = 
correspond to the solution of Schrodingers equation for the two particle problem. Note that we can 
now also determine via Eqs. (j55p and (j49p the ^-dependent matrix Q in the Hubbard-Stratonovich 
transformation ()43p . 

To summarize, we found the solution to the flow equation by first adjusting the scale-dependent 
Hubbard-Stratonovich transformation dkQ~^ such that no additional term is generated by the 
flow. All contributions to an interaction of this form are dynamically expressed by bound state 
exchange processes. As a second step we choose a A;-dependent linear field redefinition (encoded by 
the matrix M^) such that also the Yukawa interaction /i$ remains ^-independent. While the first 
step (to enforce vanishing A^) is rather natural, the second is more arbitrary. In our case the choice 
of Mk was guided by the insight we obtained in section IIIII More general it is useful to redefine 
the composite fields such that the Yukawa-like interaction /i$ remains independent of the center 
of mass momentum p, which is always possible. For some problems and for suitable choices of the 
cutoff function i?^ it may also be possible to choose the basis for the composite fields <I> such that 
the inverse propagator P<j> is a diagonal matrix. As discussed in section IIIII this is possible for the 
problem considered here for very simple cutoff functions such as Rk = k"^- For more complicated 
Rk one can still diagonalize P$ but the drawback is then that might depend on p. 



V. APPROXIMATIVE SOLUTIONS 

Up to this point we discussed the application of the functional RG to the nonrelativistic two- 
body problem in a rather formal way. No approximations were needed, but we had to introduce 
an additional functional integral over a bilocal composite fields ^{x,y) where y labels the relative 
coordinate of the two particles. For more complicated problems it is difficult to follow the RG flow 
for bilocal fields. It is therefore necessary to find useful approximation schemes. One possibility is 
to expand the field in terms of an complete orthonormal set of functions fn{y) like 

^{x,y) = Y,<t>n{x)Uy). (58) 

One can then consider the flow of the propagator and coupling constants in that basis. A sensible 
approximation would now be to take only a finite subset J of the infinite index set X into account 
and to neglect all couplings between the fields (pn, n £ J and the fields (pm n ^ J^- The approx- 
imation becomes good if the infiuence of the neglected couplings onto the flow of the considered 
quantities is small. This will, of course, only work for some particular choices of the set fn{y)- 
In the following we will consider as an example a Yukawa potential of the form 

V{X -y)= ,^ ^ -m\xr-X2\_ (59) 
47r|xi — X2I 

Note that for m — t- 00 this potential approaches a contact interaction 

V{x-y)^ ^5^^\x-y). (60) 

We choose a particular simple set of orthonormal functions 

$(x,y) = <Pr.lm{x) Yi^{ng) RnlM). (61) 

n,Lm 
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Here we use the spherical harmonics and some functions Rni for the radial direction. A 
possible choice would be the associate Laguerre Polynomials or any other suitable normalized set 
of orthogonal functions. Depending on the cutoff function, it may be possible to find analytic 
solutions for the flow equations in a particular basis. 

As a test of the robustness of our formalism, we will try a very crude approximation in the 
following and include only the term with / = m = in Eq. ()6ip . Also the radial dependence is 
truncated to a single function R{\y\) with support only for \y\ = 0. In other words, we use 



It is useful to introduce the real auxiliary field a which can be seen as the remnant of a massive 
photon in the nonrelativistic limit. Similar to the Hubbard-Stratonovich transformation performed 
in section IIIIl we multiply to the partition function a functional integral over the fields cJi and a2 
weighted by the quadratic action 




(62) 



The microscopic action corresponding to Eq. (j59p reads in momentum space 




xilJi{qi)'ip2{q2)S{qi + ^2 - 93 - Q^)- 



(63) 




(64) 



with 




g 



(65) 



and similar for {ijj2'4'2)- Adding this to Eq. (|63|) we arrive at 




(66) 



+0-2 (pi -p2)V'l(Pl)V'l(?'2)}- 



(67) 
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For the flowing action we make the fohowing truncation 



Jp 




il^limp - (iiWiimp + qi) 



(68) 



X X^{p) ipiirjip + q2)tp2{'n2P - <12)- 



We win use our freedom in the choice of the scale-dependent Hubbard-Stratonovich transformation 
to ensure that A^(p) = at all scales. It is a consequence of the nonrelativistic dispersion relation 
that the couplings in the a-sector of the theory (the charge e and the propagator of the u-boson) 
do not receive any modifications from the renormalization group flow. This can also be checked 
explicitly by looking at the flow equations for these quantities. One advantage of introducing the 
field a instead of working with V{q — (f) as before is that one can also introduce a cutoff function 
for fj. In principle, one could do this also for the composite boson (p, but we will not do this for 
simplicity, here. We use 



with Rk{z) = - z)e{k'^ - z). 

Now that we have fixed the truncation and the cutoff function, we can determine the flow 
equations. By setting the ansatz in Eq. (j68p into the flow equation (|48p we find in the center of 
mass frame {p = 0) 




+ai{-p)[-Rk{f)]a2{p) 



(69) 




1 



(70) 



po-^,(.(f+Rkm) 

+h^{pfdkQ-Hp), 




1 



(71) 



po-ii'f+Rkm) 

~p^{p)h4p)dkQ-Hp) - h^{p)dkM{p) 




P4>{P?dkQ'\p) - 2h^{p)dkM{p). 



(72) 



Here, the derivative hits only the explicit A;-dependence of the cutoff function R^. Before we 
solve the flow equations it remains to perform the integration over the spatial momentum q. For 
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the cutoff function in Eq. ()69p the integration is very simple and gives 



+h^{pfdkQ~\p), (73) 



67r2 (m^ + k^){po — j^k^) 
PMh4>{p)dkQ-\p) - h^{p)dkM{p), (74) 



OkP^ip) = — — ^fc 



-P^{p)^dkQ-\p) - 2P^{p)dkM{p). (75) 

We choose now the scale-dependent Hubbard-Stratonovich transformation dkQ~^{p) such that 
dk^iiip) = 0. This gives 

dkQ-\p) = -6^2;,^(p)2^'=(^2 + p)2(p^__Lfc2)- (76) 

In addition, we choose the rescaling of the field (p{p) encoded by dkM{p) such that dkhfp{p) = 
dkh(f,{Q). What remains then is the freedom to choose (9fcM(0). This is just the freedom to make 
a p-independent rescaling of the field (j). Usually, one chooses this wavefunction renormalization 
such that the propagator for the field has a residue of value unity at a pole that corresponds to 
a propagating particle. However, since dkM{0) can be chosen arbitrary in principle, we choose it 
here for simplicity such that 5^/1^(0) = 0. This results in 

''^'^'^ = -'^' \m^ + k-)lpo-^^k^) ^''^ 

^1^^'^^^^^^ (m2 + k^Yip, - l^k^) ■ 

We keep in mind, that with this choice the poles in the propagator for the field (j) may have residues 
with values different from one. To summarize, we have obtained the flow equations 

dkMv) = 0' (78) 
dkh4,{p) = 0, (79) 
hlk^ 



+ 2P^{p)—^dk 



1 



e 



k^ . 1 



- ^^^^)'/5^^'=(..2 + P)2;po-^it2)- (80) 

Before we proceed with the solution of the flow equation let us discuss the initial conditions 
for large values of the scale A; — )• A. In the limit — )• A — )■ cx) the contribution of the i^-particle 
exchange to the effective interaction between the '0-particles should vanish which implies 

lim = 0. (81) 

A-j-oo A 
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For large but finite A this value gets modified by one-loop contributions which can be obtained 
from the flow equation for 

_ f 



P^,A Jg {m^ + g^ + Rum? Po - + ^a(^"")) 



37r2A3 ■ 



(82) 



In the last line we assumed S> m'^, A^ S> 2/ipo- Indeed we find that this vanishes in the limit 
A — 7- oo. Since only the above ratio is fixed, there is still some ambiguity in the choice of the initial 
conditions. Indeed it is always possible to rescale the field (p which changes /i| ^ and P^^a but keeps 
the ratio fixed. Physical observables are of course not affected by such a rescaling. For simplicity 
we choose 

h(j>A = = e, 

S7r2 

P^,A = = -AVc. (83) 

The last equation defines the abbreviation c. We have now all ingredients to solve the flow equation. 
From Eq. (jSOp we obtain 

dkP. 



^ 2cP, 



(1 - 2//po/fc2)2 

1 



A;2 



+ 



1 



(1 - 2/ipo/A;2)(l +m2/A;2)2 
cPl 



(84) 



fc4 



1 



+ 



(l-2/ipoA2)2(l + m2/A;2)2 
2 



(l-2/ipoA2)(l+mVA:2)3 
We start our investigation with po = and m'^ = 0. The flow equation ([8^ simplifies to 



4cP^ 3cP^ 

dkP^ = -c+-gf--^. (85) 

The solution of this flow equation with the initial condition in Eq. (jSSp is quite simple. It is given 
for A — )■ oo by 

P^ = -k^/c + k^. (86) 

We note that the propagator for the composite field cj) changes its sign at the scale k = c. The zero 
crossing of the iinverse propagator corresponds to a bound state. Would we solve the flow equation 
for a negative frequency po < 0, would the terms involving pQ in Eq. ()84p act as an infrared cutoff 
and effectively stop the flow at the scale where fc2 « 2//|po| ■ Choosing |po| ~ c^/ (2/^) leads to P,^ = 
at the macroscopic scale A; = which corresponds to a pole in the propagator for the composite 
field 4>. We note that c^/ (2^) is proportional to the Rydberg energy ~ e^/i and the bound state can 
therefore be interpreted as the lowest bound state of the Coulomb potential. That this comes out 
qualitatively correct is encouraging, although quantitative agreement with the expectations from 
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quantum mechanics can not be expected in this case due to the simple approximation made in Eq. 

Let us now come to the limit of very large m^. The potential becomes now of the singular form 
in Eq. (|60p . For that reason, the theory must be regularized by introducing a finite UV cutoff 
Auv- However, we will see that many results do not depend on the precise choice of Auv- The 
initial values of couplings are fixed by imposing h'^\/P(j),A = at A = Auv- In addition we choose 

= e as before. The flow equation ()84p becomes for large 

= f-^ + ^ - ^1 \,2,2 (87) 
with the abbreviation w = —2fj.pQ. It has the solution 



(88) 



- 2m^{k'^ + w)i^- ck{2k'^ + 3w) 

+ {k^ + w) [Sc-v/w aTi:ctan{k / ^/w) + D] 

with an arbitrary constant D. The condition that e^/P^^A = fixes 

^ cA(2A2 + 3w) r- / . , ^^ 

D = — ^— ^ - 3cy/w arctanfA/Viy). (89) 

We note that for large k but k < A the inverse propagator is negative. Depending on the choice 
of the parameters and (or c) may change its sign during the flow towards small values 
of k. This corresponds then again to a bound state. To see this we consider P^ for /c = which 
becomes for A^ S> w 



m 



4 



P^{k = 0) = m'- — 3-^. (90) 



:(A 



4 



To interpret this expression let us assume in the following that the fields V'l and V2 describe 
fermions with equal mass Mi = M2 = M, ^ = M/2. We can obtain the scattering length a from 
Eq. dHO]) as 

M 

a = -^X^,eS, (91) 

with 

A^^.eff = 2 + ~D~n n n\ = ~\ 2 ' (^^) 

F^[k = U,w = U) cA — 



For a more detailed explanation of Eq. (|9ip we refer to refs. 0, [Sj- Insert ing the abbreviation 
c = 2^e^/(37r^) this gives 

~ 4A _ 47rm^ ' (^^) 

Depending on the choice of the initial parameters and for a given value of A, the scattering 
length can be both negative and positive. For large positive scattering length one expects the 
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presence of a shallow dimer state in the spectrum 16|]. Its binding energy can be calculated from 



the on-shell condition of the i?i)-particle which implies P,/,(/c = 0) = 0. From this we find 



A-^V^)=m2 (94) 



or after reinserting w = — 2//po = 2^£' 



This is precisely the well known relation between the scattering length a and the binding energy 
E of the shallow dimer that exists for a > [3] ■ 

It is interesting that our treatment yields exact results in the limit of large m. On the other side 
this could have been expected, since a quite similar formalism using a fc-independent bosonization 



yields exact results for the contact potential, as well [17Hl9l|. In any case it is interesting that 
the relatively simple approximation including only one boson already yields qualitatively correct 
results for the Coulomb potential and quantitatively precise results for the contact potential. More 
elaborate approximations will allow quantitative investigations for a large class of nonrelativistic 
and instantaneous interactions. Since this is not the purpose of the present paper we will discuss 
in the last chapter some generalizations of the formalism to problems that are in a certain sense 
more complicated than the nonrelativistic few-body problem. 



VI. GENERALIZED FORMALISM 

While most parts of this paper where devoted to the treatment of nonrelativistic particles with 
instantaneous interactions, we will use the present section to present some generalizations of the 
formalism to a wider class of problems. For this purpose we first somewhat generalize the notion 
of a composite field. For some set of functions fn{y) where y = {yo,y) is a space-time coordinate 
and for two real numbers r/i, r/2 G (0, 1) with rji + rj2 = 1 we define 

{ipii^2){x,y) = Tpi{x + r]2y)i^2{x - rjiy), 

ii^Mnix) = [ f*{y){^i^2){x,y). (96) 

We assume that the functions fn{y), n €z M constitute a complete orthonormal set such that the 
following relations hold 

fn{y)fm{y) = Snm, 

^fn{^)fn{y) = 6{x-y). (97) 

neJ\f 

In momentum space with fn{y) = /p e~*^^/n(p), /^(y) = /pe*^^/*(p) these relations become 

fn{p)fm{p) = Snm, 

^fn{p)fn{q) = 6{p-q), (98) 
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and the composite field reads in momentum space 

(V'lV'2)n(p) = / fni<l)'^liVlP + Q)miV2P- q), 
Jq 

{^Mlip) = [ fniq)r2{V2P-q)ri{ViP + q)- (99) 



q 



The index n we use to label the functions is an abstract index that can be either continuous 
or discrete. Note that the bilocal fields used for the investigation of nonrelativistic particles with 
instantaneous interactions can be embedded into the above description by choosing fn{y) to have 
support only for i/q = 0, i.e. 

fn{y) = S{yo) gn{y). (100) 

Similar to previous chapters we introduce bosonic fields for every composite combination. For 
example, the term in the action corresponding to the Yukawa-type coupling is 

^i^U) =-Y^j {<PUp)hip)nn.ii^li^2)M + C.C.} . (101) 

It will often be possible and convenient to work with a Yukawa coupling that is independent of p 
and diagonal with respect to the indices n, m. The action involves also a term quadratic in the 



fields (/)*, (j)n- We write it as 



= E / ^*m{P)P<t>{P)nm(t^M- (102) 
n,iri P 

Similarly, we formally introduce also a quartic interaction term for the fields ip 

j.(^,4) ^ _ ^(V;iV2):(p)A^(p)„„(VlV'2)m(p). (103) 

n,m 

As an important side remark, we note that a very large class of quartic interactions can be 
writte in the form (|lU3p . A completely general interaction between the fields ipi and ip2 can be 
written in a homogeneous situation as 

Jki...k4 (104) 

V'i(^3)V'2(A;4) ^{h + k2- h- ki). 

Due to momentum conservation, the vertex function depends actually only on three independent 
variables. One possibility to choose them is 

p = ki + k2 = ks + k^, 

qi = m^i - mk2, 

q2 = mh - mki- (105) 

Eq. (jl04p can then be written as 

ri'^''^ = -/ (V'lV^2)*(p,(Zl) A^(p,(?i,(?2) (VlV'2)(p,g2), (106) 
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with 



(107) 



Assuming now that X^{p, (71, 92) is regular with respect to the arguments qi and (72 we can expand 

) = Yl ^Ap)nmfmiQ2), (108) 



with coefficients 



Av.(p)r 



fnin) A^(p, 91,^2) fm{q2)- 



(109) 



91.92 



Plugging Eq. (jlOSp into Eq. ()106p and using ()99p we find that this is precisely of the form in Eq. 
()103p . Note that we did not assume that A^(p, 51,(72) is regular with respect to the argument 
p. Possible singularities are transfered to the coefficients \^{p)nm- In fact, it is a big advan- 
tage of our formalism that poles or branch cuts in \^{p)nm can be described by relatively simple 
parameterizations of the inverse propagator P^{p)nm for a composite particle. 

Let us assume that for a fixed Hubbard-Stratonovich transformation a term of the form in Eq. 
(|103p is generated by the renormalization group flow, i.e. 9fcA^(p)„m| jjs 7^ ^^"^ then follow 

the calculation in Sect. IIVI and use a scale-dependent Hubbard-Stratonovich transformation to 
realize dk\^{p)nm = 0. This will lead to additional contributions in the flow equations for h^{p)nm 
and P(f,{p)nm- In addition, a fc-dependent linear transformation of the fields (/>n(p) can be employed, 
for example to keep h^{p)nm independent of p. More concrete, the flow equations can be derived 
similar to Eq. (|52p and read in symbolic notation 

+ h^{p)nrdkQ~^ip)rsh<f,{p)sm, 

r,s 

dkh(l,{p)nm = 5fc/l<^(p)„m|jjg 

-y^,P'i>{p)nrdkQ^^{p)rsh^{p)sm 
r,s 

-'^idkMl) (p) nr (P) rm , 
r 

dkP<j,{p)nm = dkP^{p)nm\^g 

- ^ P4>{p)nrdkQ~^{p)rsP4>{p) rm 

r,s 

- Y{dkMl){p)nrP<f>{p)rm 

r 

-Y,PMnr{dkMk){p)rm. (HO) 

r 

We see that by choosing dkQ~'^{p)nm conveniently, we can absorb a quartic interaction term A^ 
as in Eq. ()104p that is regular with respect to the relative momenta qi and q2 into the exchange 
of a composite boson of the type ipiip2- Besides this particle-particle pair also other combina- 
tions are possible, such as the combination ViV'i or V'2V'2 corresponding to particle-hole pairs in 
nonrelativistic physics or a particle-antiparticle pair in relativistic quantum fleld theory. With a 
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construction similar to the one presented above, one can show that a contribution to the general 
interaction term in Eq. ()104p that is regular as a function of ki + kj, and k2 + ki^ (but might 
have singularities with respect to ki — kj, = ki^ — ^2) can be absorbed into the exchange of a V'lV'i 
pair. The construction would be similar to introducing the cr-field in Eq. ()64p but now with a 
A;-dependent Hubbard-Stratonovich transformation. In a relativistic field theory (or in a nonrela- 
tivistic field theory at nonzero density) the propagator of this field might get a nontrivial frequency 
dependence and thus become a dynamical field. Finally, the exchange of a '(/'i'02-pair can describe 
contributions that are regular as functions of the combinations ki + ^4 and k2 + k^ but might have 
singularities with respect to ki — k^ = k^ — k2. 



VII. CONCLUSIONS 

We have shown in this paper how a recently derived exact flow equation can be employed for the 
treatment of bound state formation in quantum field theory. For a nonrelativistic field theory, the 
relevant flow equations can be integrated exactly and our approach is equivalent to the standard 
quantum mechanical treatment. The presented formalism can also be employed to find approximate 
solutions to bound state problems. This can be useful under circumstances where no exact solution 
from quantum mechanics is available. Examples for such situations arise in relativistic quantum 
field theory, due to non-instantaneous interactions or in nonrelativistic quantum field theory at 
nonzero density and temperature. 

One big advantage of our formalism is that fundamental and composite particles are described 
in a very similar way by fluctuating quantum fields. A scale dependent Hubbard-Stratonovich 
transformation allows to absorb complicated interaction vertices into relatively simple tree diagrams 
of a Yukawa-type theory. Depending on the situation, a particular composite field might become 
a propagating degree of freedom, end as a gapped excitation or remain an auxiliary field without 
own dynamics when the infrared cutoff scale is lowered. The formalism is particularly useful in 
situations where the relevant degrees of freedom change during the renormalization group fiow. 
This is often the case when interaction effects are strong. Prominent examples are QCD where 
the ultraviolet physics is dominated by quarks and gluons while mesons and baryons dominate the 
physics in the infrared or the Hubbard model, where fermionic degrees of freedom dominate the 
UV, while various competing pairing instabilities dominate the infrared physics. For the example 
of QCD the fiow from UV to IR degrees of freedom was discussed in a setup closely related to the 
one discussed here 0, 21]. 



Functional renormalization group equations have been applied with quite some success to many 
problems in quantum field theory. For reviews see r ef. [lM3l. It is a general feature that 
attempts to increase the precision of these calculations face the difficult problem to find an efficient 
parameterization for the momentum dependence of vertex functions. The formalism presented in 
this paper does not only provide a general framework for such parameterizations (similar to the 
one proposed in ref. [22]) but also allows to take the essential composite degrees of freedom as 
dynamical fields into account. This includes also possible interactions between composite particles 
which typically correspond to complicated and nonlocal higher-order interactions in terms of the 
fundamental fields. Another advantage of the description involving the composite fields is that 
spontaneous symmetry breaking, for example in a superconducting phase, is straight-forward to 
describe. 

Many of the phenomena discussed above - including the formation of bound states - can also 
be described with other functional methods of QFT such as Schwinger-Dyson and Bethe-Salpeter 
equations. Renormalization group equations have the advantage that they allow a straightforward 
discussion of interesting features such as fixed point behavior and universality. Close to renor- 
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malization group fixed points one can distinguish between relevant and irrelevant (and marginal) 
operators. This implies that the flow equations determine which terms in a given expansion become 
important and which can be neglected! Often, a large class of microscopic theories are attracted 
towards the same fixed point and the system shows universal features that are independent of the 
concrete microscopical realization. Using a scale-dependent Hubbard-Stratonovich transformation, 
one can study for a theory with bound states the different fixed points as well as the crossovers 
between them |20|]. 

In summary, we have developed a formalism to describe bound states in quantum field theory 
using an exact flow equation and look forward to applications to many interesting problems in the 
future. 
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